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The microscopic description of nuclear fission still remains a topic of intense basic research. Un¬ 
derstanding nuclear fission, apart from a theoretical point of view, is of practical importance for 
energy production and the transmutation of nuclear waste. In nuclear astrophysics, fission sets the 
upper limit to the nucleosynthesis of heavy elements via the r-process. In this work we initiated 
a systematic study of intermediate energy proton-induced fission using the Constrained Molecu¬ 
lar Dynamics (CoMD) code. The CoMD code implements an effective interaction with a nuclear 
matter compressibility of K=200 (soft EOS) with several forms of the density dependence of the 
nucleon-nucleon symmetry potential. Moreover, a constraint is imposed in the phase-space occu¬ 
pation for each nucleon restoring the Pauli principle at each time step of the collision. A proper 
choice of the surface parameter of the effective interaction has been made to describe fission. In 
this work, we present results of fission calculations for proton-induced reactions on : a) 232 Th at 27 
and 63 MeV, b) 235 U at 10, 30, 60 and 100 MeV, and c) 238 U at 100 and 660 MeV. The calculated 
observables include fission-fragment mass distributions, total fission energies, neutron multiplicities 
and fission times. These observables are compared to available experimental data. We show that 
the microscopic CoMD code is able to describe the complicated many-body dynamics of the fission 
process at intermediate and high energy and give a reasonable estimate of the fission time scale. 
Sensitivity of the results to the density dependence of the nucleon symmetry potential (and, thus, 
the nuclear symmetry energy) is found. Further improvements of the code are necessary to achieve 
a satisfactory description of low energy fission in which shell effects play a dominant role. 

PACS numbers: 25.70.-z, 25.70.Hi,25.70.Lm 

Keywords: proton-induced fission, fission fragments, mass yield distributions, neutron multiplicities, fission 
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I. INTRODUCTION 


Nuclear fission, since its discovery more than 75 years 
ago, has played a prominent role in applications as well 
as basic nuclear research. Among the wide range of mod¬ 
ern applications and given the increasing energy demand 
worldwide 1L nuclear energy production in Generation- 
IV reactors 2, 3} and the incineration of nuclear waste 
in Accelerator-Driven Systems (ADS) 0,0 are areas of 
current intense efforts. In parallel, fission offers an im¬ 
portant mechanism to produce a variety of isotopes for 
medical and industrial use {6:j. It is also one of the main 
approaches to produce exotic neutron-rich nuclei in Rare 
Isotope Beam (RIB) facilities 0®. 

From an astrophysical point of view, fission is a key 
reaction of the rapid neutron capture process (r-process) 
and essentially sets the u ppe r boundary on the synthe¬ 
sis of heavy elements jl3l - ll71 . In the same vein, fission 
largely determines the stability and the properties of su¬ 
perheavy elements (l8l - l2l| . 


* Corresponding author. Email: soulioti@chem.uoa.gr 


Understanding the mechanism of nuclear fission, that 
is the transformation of a single heavy nucleus into two 
receeding fragments, has been a long journey of fruitful 
research and debate and, still today, is far from being 
complete. Upon its discovery, fission was interpreted by 
Meitner and Frisch 22] as the division of a charged liquid 
drop due to the interplay of the repulsive Coulomb force 
between the protons and the surface tension due to the 
attractive nucleon-nucleon interaction. In the seminal 
paper of Bohr and Wheeler [23}, fission was described 
with a liquid-drop model and the first estimates of fis¬ 
sion probalilities were obtained based on statistical argu¬ 
ments. The first detailed calculations of potential ener¬ 
gies of deformed nuclear drops were performed by Frankel 
and Metropolis in 1947 [24| employing the ENIAC, one 
of the first digital computers. Despite the success in the 
interpretation of fission based on the liquid-drop model, 
the prevailing asymmetry in the mass distribution of the 
minor actinides could not be deciphered until shell cor¬ 
rections to the macroscopic liquid drop were taken into 
account (see below). 


A detailed statistical model that could describe asym¬ 
metric fission was developed by Fong [25|. Further ad¬ 
vancements lead to the scission model of Wilkins [26| and 
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the random neck-rupture model of Brosa f27t . A current 
version of the latter is the temperature-dependent Brosa 
model developed in [28}. These models offer main ingredi¬ 
ents in current statistical models of fission (e.g. [29H-!.'Sl | 1. 
Two widely used current models that contain a statisti¬ 
cal description of fission are the code GEMINI++ [34] 
in which fission is described by the transition-state ap¬ 
proach of Moretto [35[ and the code SMM (j36j, and refer¬ 
ences therein) in which low-energy fission is described by 
an empirical parametrization and higher-energy fission 
is treated as one of the possible channels of statistical 
multifragmentation. 

Along with the statistical description of fission, the 
dynamical approach to the fission process was put for¬ 
ward early on in the seminal paper of Kramers [37} . 
Based on the assumption that the deformation degree 
of freedom can be viewed as a Brownian particle inter¬ 
acting stochastically with the single-particle degrees of 
freedom that constitute a heat bath, Kramers analyti¬ 
cally solved the appropriate ID Fokker-Planck equation 
and predicted that fission is actually retarded relative 
to its rate obtained from purely phase-space arguments. 
Along these lines, dynamical approaches based on the 
dissipative character of the nuclear shape motion were 
developed based either on the Fokker-Planck equation 
or the Langevin equations in more that one demensions 
describing the deformation degrees of freedom. A de¬ 
tailed review of these efforts is given in [H, [HJ. Refined 
dynamical approaches along these lines have continued 
with increasing degree of sophistication and success (e.g. 

EM3). 

Notable to all efforts to descibe fission based on the 
macroscopic LDM is the prediction of symmetric mass 
yields, in contrast to the large body of experimental data 
of minor actinides that indicate asymmetric mass yield 
distributions of low-energy fission. This discrepancy was 
remedied by Swiatecki [571 ] by the inclusion of a micro¬ 
scopic correction to the macroscopic LDM part of the po¬ 
tential energy surface (PES) of the nucleus. A detailed 
shape-dependent macroscopic-microscopic PES was ob¬ 
tained by Strutinsky [58l, [HI] forming the basis of the 
successfull shell-correction approach to the PES. 

The Strutinsky macroscopic-microscopic approach, 
due to its simple physical foundation and numerical flexi¬ 
bility, has seen continuous development and success. De¬ 
tailed relevant reviews concerning the description of fis¬ 
sion barriers can be found in [60H62| . The most detailed 
prediction of PES based on the macroscopic-microscopic 
approach was performed by Moller et al. on a five- 
dimensional deformation space [60, [63]. This detailed 
PES description forms the basis of the recent dynami¬ 
cal description of fission in the limit of strong dissipative 
coupling (Smolutsovski limit) in which the fission process 
resembles a random walk on the nmltidimentional PES 



Apart from the macroscopic-microscopic description of 
the heavy-element PES relevant to the description of nu¬ 
clear fission, fully microscopic approaches based on the 


nuclear Density Functional Theory (DFT) have been de¬ 
veloped. Some recent representative works are 
in which extensive reference to previous works can be 
found. We also report the recent work [73| on uncon¬ 
strained DFT calculations in which properties of very 
deformed nuclei pertinent to fission are described. 

Whereas the static properties of the PES are very well 
accounted for by the modern DFT-type approaches, the 
microscopic description of the full dynamics of the fis¬ 
sion process still remains a daunting project for nuclear 
theory. 

Two main quantal approaches have been adopted in 
the past to describe the fission dynamics. First in 1978, 
the time-dependent Hartree-Fock (TDHF) theory was 
applied to the fission process Q. However, as docu¬ 
mented in recent studies of heavy-ion collisions [zlHzJ, 
the TDHF approach, being a one-body (mean-field) ap¬ 
proach has essentially no correlations beyond the mean 
field and, as such, cannot fully describe fluctuations (as, 
e.g., encountered in nucleon exchange or fragment for¬ 
mation). The lack of fluctuations also affects the correct 
triggering of scission in the time evolution of a deforming 
nucleus. 


The second approach is based on the adiabatic approx¬ 
imation and involves the adiabatic TDHF theory [78} and 
the time-dependent generator coordinate method (GCM) 
7!H . In these approaches, the adiabadic hypothesis is in¬ 
voked for the fission process, namely, that the nucleonic 
degrees of freedom are fully equilibrated during the slow 
evolution of the deforming nucleus over the macroscopic 
PES determined by the appropriate deformation degrees 
of freedom. We note that the adiabatic approximation is 
also inherent to the Langevin (or Fokker-Plank) type of 
approaches mentioned earlier. Under this approximation, 
the dynamics is described up to the scission configura¬ 
tion. At this point a sudden approximation is invoked (as 
in the early scission model [2g]) in order to obtain the fis¬ 
sion fragments and their characteristics (mass, charge, ki¬ 
netic energy). However, as stated in recent works [80118l| . 
ignoring the non-adiabatic effects during and past scis¬ 
sion leads to limitations in the predictive power of the 
models regarding the fragment characteristics (e.g., the 
kinetic energy). 

The recent dynamical study 81} attempts to bridge 
the two regimes. The adiabatic phase of the fission pro¬ 
cess is described with a static mean-field (DFT) approach 
and the nonadiabatic phase is carefully described with 
TDHF. It is found that the proper treatment of the nona¬ 
diabatic phase results in an accurate description of the ki¬ 
netic energy and excitation energy of the resulting fission 
fragments. Before closing the above review on dynami¬ 
cal approaches, we also report efforts to describe sponta¬ 
neous fission employing TDHF and Feynman’s path inte¬ 
gral approach [82} , as well as the semiclassical equivalent 
method employing Vlasov’s equation 


The preceeding rather limited overview of the extended 
literature on nuclear fission dynamics clearly indicates 
that a full microscopic description of the fission pro- 
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cess is unusually challenging and, as of today, “it has 
not been possible to establish a computationally feasible 
framework capable to describe real nuclei with realistic 
interactions” j64|. Nevertheless, we can see that sub¬ 
stantial progress in both macroscopic, as well as quantal 
dynamical approaches continues vividly. Motivated by 
the current situation regarding fission dynamics, in the 
present work, we initiated a study of fission based on the 
semiclassical microscopic N-body Constrained Molecular 
Dynamics (CoMD) model [85| in regards to its ability 
to describe the full dynamics of the fission process in 
proton-induced reactions on Th and U from low to high 
energies. 

In the following discussion, we classify, as customary, 
the proton-induced fission reaction according to the pro¬ 
ton energy E p as: a) low-energy when E p < 20 MeV, 
b) intermediate energy when 20 MeV < E p < 200 MeV 
and c) high energy when E p > 200 MeV. The present 
paper has the following structure. In section II, we high¬ 
light the basic aspects of the CoMD code and present 
the way that the code is applied to nuclear fission. In 
section III, we present results of fission calculations for 
proton-induced reactions on : a) 232 Th at 27 MeV and 
63 MeV, b) 235 U at 10 MeV, 30 MeV, 60 MeV and 100 
MeV, and c) 238 U at 100 MeV and 660 MeV. The fis¬ 
sion fragment mass yield distributions are presented, as 
well as total fission cross sections, total energies, neutron 
multiplicities and fission times with respect to incident 
proton energy. In section IV, discussion of the results 
and conclusing remarks are given. 


II. DESCRIPTION OF THE THEORETICAL 
MODEL 

The theoretical model employed in this work is the 
microscopic Constrained Molecular Dynamics (CoMD) 
model originally designed for reactions near and below 
the Fermi energy [85l - l87| . Following the general approach 
of Quantum Molecular Dynamics (QMD) models [88j, in 
the CoMD code nucleons are described as localized Gaus¬ 
sian wave packets. The wave function of the N-body 
nuclear system is assumed to be the product of these 
single-particle wave functions. With the Gaussian de¬ 
scription, the N-body time-dependent Schodinger equa¬ 
tion leads to (classical) Hamilton’s equations of motion 
for the centroids of the nucleon wavepackets. The po¬ 
tential part of the Hamiltonian consists of a Skyrme-like 
effective interaction and a surface term. The isoscalar 
part of the effective interaction corresponds to a nuclear 
matter compressibility of K=200 (soft EOS) or K=380 
(stiff EOS). For the isovector part, several forms of the 
density dependence of the nucleon-nucleon symmetry po¬ 
tential are implemented. Two of them will be used in the 
present work: the “standard” potential [red (solid) lines] 
and the “soft” potential [blue (dotted) lines] in the figures 
that follow. These forms correspond to a dependence of 
the symmetry potential on the 1 and the 1/2 power of 


the density, respectively (see, also, [89} and references 
therein). 

We note that in the CoMD model, while not explicitly 
implementing antisymmetrization of the N-body wave- 
function, a constraint in the phase space occupation for 
each nucleon is imposed, effectively restoring the Pauli 
principle at each time step of the (classical) evolution of 
the system. This constraint restores the fermionic na¬ 
ture of the nucleon motion in the evolving nuclear sys¬ 
tem. More specifically, at each time step, and for each 
nucleon, the presense of neighboring nucleons is deter¬ 
mined in phase space. If the phase space occupation 
probability is greater than 1, then the code changes the 
direction (not the magnitude) of the nucleon momentum, 
so that the total momentum and kinetic energy are con¬ 
served fH}. 

The short range (repulsive) nucleon-nucleon interac¬ 
tions are described as individual nucleon-nucleon colli¬ 
sions governed by the nucleon-nucleon scattering cross 
section, the available phase space and the Pauli princi¬ 
ple, as usually implemented in transport codes (see, e.g. 
[901]'). The present CoMD version fully preserves the to¬ 
tal angular momentum (along with linear momentum and 
energy), features which are critical for the accurate de¬ 
scription of observables from heavy-ion collisions and, for 
the present study, the fission dynamics. 

The ground state configurations of the target nu¬ 
clei were obtained with a simulated annealing approach 
and were tested for stability for long times (2000-3000 
frn/c). These configurations were used in the subsequent 
particle-induced fission simulations. 

In the calculations of the present work, the CoMD 
code was used mainly with its standard parameters. 
The soft density-dependent isoscalar potential was cho¬ 
sen (K=200), as well as the standard and soft symmetry 
potentials, as mentioned above. The surface term was 
set to zero. 

For a given p-induced reaction, a total of 3000-5000 
events were collected. For each event, the impact param¬ 
eter of the collision was chosen in the range b = 0-6 fm, 
following a triangular distribution. Each event was fol¬ 
lowed in time up to 15000 fm/c ( 5.0xlO” 20 sec ) The 
phase space coordinates were registered every 50 fm/c. 
At each time step, fragments were recognized with the 
minimum spanning tree method [§5[ , and their properties 
were reported. From this information, we obtained infor¬ 
mation on the evolution of the fissioning system and the 
properties of the resulting fission fragments. We mention 
that we consider as fission time (t fission) the time inter¬ 
val between the implantation of the proton in the target 
nucleus and the emergence of the two fission fragments. 
We allowed an additional time of 2000 fm/c after scission 
for the nascent fission fragments to de-excite. (We varied 
this time interval from 2000 to 5000 fm/c and we did not 
notice an appreciable change in the characteristics of the 
fission fragments.) Thus, in the following discussion, for 
each event, the fission fragment properties are reported 
and studied 2000 fm/c after scission. 
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A typical time evolution of a fissioning system as pre¬ 
dicted by CoMD is presented in Fig. 1. The figure refers 
to p-induced fission of 232 Th at 63 MeV and gives a three 
dimensional representation of the fissioning system in the 
center of mass at three time instants. At 0 frn/c (Fig. 
la), the proton approaches the target nucleus 232 Th. At 
2000 fm/c (Fig. lb) the nucleus has been substantially 
deformed. This configuration is near or past the saddle 
point. (For this event scission occurs at t/i SS j O n=2500 
fm/c.) At 4000 fm/c (Fig. lc), we observe the two fis¬ 
sion fragments and the emission of two neutrons depart¬ 
ing nearly perpendicular to the fission axis. 

In Fig. 2, the time evolution of the axial quadrupole 
moment Q 20 of the fissioning system (Fig. 2a) and 
its mean radius (Fig. 2b) are presented for the same 
event. Both quantities increase with time and indicate 
the course of the system toward fission. (We note that 
the decreasing value of Q 20 after scission is attributed to 
the continuous rotation of the deformed nucleus and the 
resulting fission fragments.) 


III. RESULTS AND COMPARISONS 

A. Mass yields: Low and Intermediate Energy 

We begin our study of the behavior of the CoMD code 
with comparisons to the recent experimental data de¬ 
scribed in [29i|9lJ. First we show the CoMD calculations 
for the proton induced fission of 232 Th at energies 27 MeV 
and 63 MeV, using the standard and the soft symmetry 
potential. In Fig. 3, the mass yield distributions are il¬ 
lustrated for the reaction at 27 MeV. In the experimental 
data (full points), we observe the asymmetric nature of 
the fission mass yield, as expected for low-energy fission 
of minor actinides. We compare the data with our calcu¬ 
lations (open points with statistical errorbars) with the 
standard symmetry potential (Fig. 3a), as well as the 
soft symmetry potential (Fig. 3b). With both selections 
of the symmetry potential, we observe a merely symmet¬ 
ric mass distribution with a rather flat top. No clean 
sign of an asymmertic mass yield distribution is seen in 
the CoMD calculations. A hint for asymmetric distribu¬ 
tion may be implied in Fig. 3b. At this point, we note 
the finding of Nadtochy et al. [92j that in dynamical 
Langevin calculations of fission, dominance of asymmet¬ 
ric mass splits relative to the (expected) symmetric mass 
split may occur as a result of the dissipative dynami¬ 
cal behavior of the system. A similar suggestion that 
asymmetric fission may result from the hydrodynamical 
behavior of the system was first reported in [93j | before 
strong shell effects were considered responsible for the 
asymmetric mass yield distribution of actinides [58|. 

The main reason for the symmertic mass yield distri¬ 
bution obtained by the CoMD code is that the nucleon- 
nucleon interaction in the model does not include spin 
dependence, thus the resulting mean-field potential does 
not contain a spin-orbit contribution. Thus, the model 


does not predict the correct shell effects in the single¬ 
particle motion of the deforming nucleus, which are nec¬ 
essary to lead to the asymmetric fission of 232 Th. A closer 
inspection of the two yields calculations, shows that they 
are slightly different but neither of them tends to resem¬ 
ble the experimental distribution. 

We wish to comment that, while in the present imple¬ 
mentation of the CoMD model the interaction has no spin 
dependence, (thus CoMD cannot describe the correct 
shell effects), the code emulates the quantum behavior of 
the deforming nuclear system, thus we should expect shell 
effects (at the mean-field level) which would correspond 
to those obtained by a deforming harmonic oscillator or 
Woods-Saxon potenial (without a spin-orbit term) jMj. 
We also mention that application of the CoMD approach 
to light atoms has successfully reproduced the electronic 
binding e nerg ies, as well as electron radii revealing shell 
structure [95[. A study of shell effects in the present im¬ 
plementation of CoMD applied to nuclei has not been 
performed to date. Such a study with the present CoMD 
code and a possible extension of it with spin dependence 
(in the spirit of recent work on BUU (96[) will be under¬ 
taken by us in the near future. 

In Fig. 4, the mass yield distribution for the same 
reaction at proton energy 63 MeV is presented. It is ev¬ 
ident that the structure of the experimental mass yield 
curve tends to become more symmetric at this higher 
energy. This is to be expected, because as the proton 
energy and, thus, the excitation energy of the fissioning 
system increases, shell effects will begin to fade (see e.g. 
[Hi HI). However, it seems that this beam energy is 
not high enough to completely wash out the shell effects, 
as two asymmertic shoulders appear in the experimental 
mass yield curve. In Fig. 4a, we show the CoMD calcu¬ 
lations with the standard symmetry potential and in Fig. 
4b, the soft symmetry potential. As in the lower energy 
case (Fig. 3), the two choices of the symmetry potential 
do not lead to substantial differences on the mass yield 
shape. 

We note that at this higher energy, as the asymmetric 
mass split is attenuated and the symmetric contribution 
is enhanced, an overall improvement in the agreement be¬ 
tween our CoMD calculations and the experimental data 
is obtained. Thus, with the current implementation of 
the CoMD code, intermediate-energy fission mass yields 
may be correctly descibed. We will explore this behavior 
with other fissioning systems in the following. 

We continue our comparisons with the recent work [98j 
on the proton-induced fission of 235 U at proton beam en¬ 
ergies of 10 and 30 MeV. In Fig. 5, we present the calcu¬ 
lated mass yield distributions (open points) of the proton 
induced fission of 235 U at 10 MeV, using the two forms 
of the symmetry potential, the standard (Fig. 5a) and 
the soft (Fig. 5b). The experimental data of j98| are 
presented by closed points. We note that the yield data 
in [98) are in arbitrary units, so we multiplied the yields 
with a factor of 7, for both reactions to make them com¬ 
parable with our calculated cross section results. In the 
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experimental data, the prevalence of the asymmetric fis¬ 
sion mode is obvious for this low-energy fission reaction. 
Our calculations with the standard symmetry potential 
indicate a rather symmetric distribution (Fig. 5a). In¬ 
terestingly, the calculations with the soft symmetry po¬ 
tential indicate a hint of asymmetric fission (Fig. 5b). 
However, as already discussed above, in the absence of 
the correct shell effects from the CoMD potential, we do 
not expect to obtain the correct asymmetric mass distri¬ 
bution. In the comparisons of Fig. 5, we observe that the 
CoMD calculations result in a wider mass yield distribu¬ 
tion as compared with the data. Apart from a possible 
calculational aspect, this may also point to a limitation 
in the data toward asymmetric mass splits. 

In Fig. 6, we present the mass yield distribution from 
the p-induced fission of 235 U at proton energy 30 MeV. 
The experimental data come from [98). In the CoMD 
calculations, we again used the two forms of the symme¬ 
try potential, the standard (Fig. 6a) and the soft (Fig. 
6b). Conclusions similar to those from Fig. 5 can be 
drawn. However, going from 10 MeV to 30 MeV proton 
energy, in the experimental mass yield distributions we 
see that the peak-to-valley ratio is reduced as the sym¬ 
metric contribution increases relative to the asymmetric 
contribution. Thus, in Fig. 6, the CoMD calculations 
(yielding a symmetric mass distribution) are obviously 
in better agreement with the experimental data. 

Furthermore we explored the behavior of CoMD at 
higher energy for the same system. In Fig. 7 we show 
the CoMD calculations for proton induced fission of 235 U 
at 60 MeV, again with the two forms of the symmetry 
potential, standard (Fig. 7a) and soft (Fig. 7b). We 
compared our calculations with available experimental 
data of proton induced fission of 238 U (not of 235 U) at 
this beam energy taken from [28^ (and normalized to our 
calculated cross sections, as in Figs. 5 and 6). We note 
that the small difference in the number of neutrons of the 
fissioning systems at this higher energy is not expected 
to substantially affect the mass yield comparisons of Fig. 
7. The experimental distribution for U indicates that the 
asymmetric fission mode still prevails. This yield curve 
presents a plateau, in contrast to the Th distribution at 
comparable energy that has a symmetric peak and two 
shoulders at asymmetric mass splits (Fig. 4). Our cal¬ 
culations, as expected, indicate a symmetric peak, as we 
saw in the lower energy cases (Fig. 5, 6) of the proton 
induced U fission. 


B. Mass yields: High Energy 


We continue our investigations with the application 
of CoMD to high energy fission reactions with protons. 
We remind that the CoMD code has been successfully 
applied to the description of a large variety of nuclear 
reactions (e.g., HMqI). As a fully dynamical code, 
we expect that it may perform well also with spallation- 
type reactions with high-energy protons. 


Toward this direction, we performed calculations for 
the proton induced fission of 238 U at 660 MeV proton 
energy. The motivation comes from the importance of 
this energy range in ADS-type applications SI and, 
specifically, the r ecent exp erimental data for this reac¬ 
tion reported in [l03l - ll06| obtained by off-line gamma- 
ray techniques. In Fig. 8, we show the experimental data 
(full symbols) that despite the large experimental fluctu¬ 
ations indicate a prevailing symmetric fission mode. (We 
point out that the yield axis is logarithmic in this fig¬ 
ure). Our CoMD calculations, with the two forms of the 
symmetry potential (Fig. 8a standard, Fig. 8b soft) are 
in reasonable agreement with the experimental data near 
symmetric mass splits. 

We note that the shape of the experimental mass distri¬ 
bution is characterized by two low-yield very asymmetric 
fission components ( “super-asymmetric” fission tl06l | [). 


Our calculations show an overall symmetric curve that 
is wide enough to contain these superasymmetric mass 
splits, predicted, however, with larger cross sections than 
the data. Our calculations resemble the wide symmet¬ 
ric mass yields observed resently in high-quality mass 
spectrometric data obtai ned in inverse kinematics at 
r elativist ic energies mm . Comparing the data of 
1104 1 with the higher-energy inverse-kinematics data 


(e.g. [1071 ]). we speculate that the former data may be 
incomplete due to the nature of the measurements and 
we suggest that measurements of this very important re¬ 
action at ^600 MeV be performed in inverse kinematics 
in the same fashion as the higher-energy data. 

We wish to point out that the lower-energy fission data 
discussed above (Figs. 3-7) have been acquired with stan¬ 
dard fission on-line counters, which cannot provide in¬ 
formation on the atomic number Z of the fission frag¬ 
ments. We remind that to obtain Z information, either 
mass spectrometric tecniques (mainly in inverse kinemat¬ 
ics) or off-line gamma-ray methods have to be used. As 
already mentioned, the data of fl03l . lT04j were obtained 
with gamma-ray techniques and thus can provide infor¬ 
mation on the Z-A correlation of the observed fission frag¬ 
ments. 


In Fig. 9, we first present the experimentally observed 
mean Z (Fig. 9a) as well as the standard deviation of 
the Z distribution (Fig. 9b.) with respect to the mass 
number A of the fission fragments. Our CoMD calcu¬ 
lations [solid (red) line] show that the fission fragments 
are more neutron-rich relative to the experimental data. 
Due to the fact that, for this high-energy reaction, the 
code causes one or two protons to be emitted before 
the scission point, we made a selection concerning the 
charge of the fissioning nucleus (Z=93) so that it cor¬ 
responds to no pre-scission proton emission at the time 
of scission. With this selection, the CoMD calculations 
[dashed (blue) line] are in better agreement with the data, 
especially for the heavier fragments. From the above 
comparison for this high-energy reaction, we may con¬ 
clude that the fission fragments, as obtained 2000 fm/c 
after scission, may still contain enough excitation en- 
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ergy to further evaporate (predominantly) neutrons and, 
thus, move closer to the data in the Z-A plot. An ex¬ 
plicit de-excitation of these fragments with a standard 
de-exciation code (e.g. [H, HH) was not performed in 
the present exploratory work, but will be performed by 
us in the near future within our plan of detailed studies 
of high-energy p-induced fission of actinides. In Fig. 9b, 
our calculations of the standard deviation are higher in 
comparison with the experimental data. However, when 
the selection of the fissioning nucleus is made, so that it 
corresponds to no pre-scission proton emission (Z=93), 
the calculations are in better agreement with the data, 
despite the large fluctuations due to the limited statistics 
of the calculations after the imposed selection. 


C. Fission Cross Sections 

After the presentation of the mass yield distribution, 
which is one of the most characteristic observables of fis¬ 
sion reactions, we continue our investigation with several 
other fission observables, starting with the total fission 
cross sections. In Fig. 10, calculated total fission cross 
sections are presented as follows. The full (red) triangles 
connected with a full line refer to the p-induced fission 
of 232 Th at the two energies 27 and 63 MeV. The full 
(red) circles connected with a full line represent the p- 
induced fission of 235 U at 4 energies (10, 30, 60 and 100 
MeV). The full (red) squares connected with the full line 
show the p-induced fission of 238 U at two proton energies 
(100 and 660 MeV). The CoMD calculations correspond 
to the standard symmetry potential. We note that the 
points at E p =660 MeV are displayed at E p =160 MeV in 
Figs. 10-15, for convenience of presentation. 

The available experimental data for 232 Th are shown 
with closed (black) triangles connected with dotted lines. 
The experimental point for 235 U at E p =60 MeV is shown 
with a closed (black) circle and for 238 U at E p =660 MeV 
with closed (black) square. 

Concerning the fission of thorium, we observe that in¬ 
creasing the proton energy (and thus the excitation en¬ 
ergy of the fissioning nucleus), there is only a slight in¬ 
crease in the calculated cross section. However, the ex¬ 
perimental data show an increase of approximately 30%, 
which our calculations do not reproduce. For the fis¬ 
sion of 235 U, in the calculations we observe a jumb of 
the cross section, with increasing proton energy from 30 
MeV to 60 MeV. The experimental point at 60 MeV is 
20% larger than our calculated point. At higher energies, 
for the proton induced fission of 238 U , the total fission 
cross section is rather constant and in rough agreement 
with the experimenal data within errorbars. 

In Fig. 11, we present the ratio of the fission cross sec¬ 
tion to the heavy-residue cross section as a function of the 
proton energy. This ratio is a very sensitive observable 
for the relative importance of fission as a deexcitation 
path for the nuclei examined. 

The CoMD calculations (with the standard symme¬ 


try potential) are shown with the closed (red) symbols 
connected with solid (red) lines with exactly the same 
correspondence as in Fig. 10. In addition, CoMD calcu¬ 
lations with the soft symmetry potential are also shown 
in this figure with closed (blue) symbols connected with 
dotted (blue) lines (such calculations were not shown in 
Fig. 10 because they would nearly overlap with the ones 
shown). 

For the fission of thorium we observe an increase in the 
ratio from 27 to 63 MeV and a rather weak sensitivity 
to the choice of the symmetry potential, the ratio being 
slightly larger at the energy of 63 MeV with the choice of 
the soft symmetry potential. For the fission of 235 U, the 
ratio increases from 10 to 30 MeV and then diminishes at 
the higher energies of 60 and 100 MeV. Decerasing trend 
also exhibits the ratio for 238 U from 100 to 660 MeV, 
possibly pointing to the increasing role of fast evapora¬ 
tion processes for the more excited nuclei involved in the 
higher energy reactions. 

Finally, focusing our attention to the behavior of the 
calculated ratio with the soft symmetry potential for the 
235 U and 238 U isotopes, interestingly we observe that this 
ratio is substantially larger than the corresponding ratio 
calculated with the standard symmetry potential. Fur¬ 
ther detailed investigation of the features of this bevavior 
is in line. However, from the present work we conclude 
that the ratio is a rather sensitive observable of the den¬ 
sity dependence of the nucleon-nucleon symmetry poten¬ 
tial, and thus, of the nuclear symmetry energy, which is 
a topic of current importance in reg ards to studies of the 
nuclear equation of state (e.g. [nmnl). 


D. Total Fission Kinetic Energy 

In the following we will examine the mean total kinetic 
energy of the fission fragments as a function of the pro¬ 
ton energy for the studied reactions. This is an important 
kinematical observable characterizing on average the de¬ 
gree of deformation, the compactness and the asymmetry 
at scission of the fissioning system and offers an impor¬ 
tant testing ground of the overall dynamical description 
offered by the employed code. In Fig. 12, we illustrate 
the mean total kinetic energy of the fission fragments for 
the aforementioned fission reactions. The symbols corre¬ 
spond to the same reactions as in Figs. 10 and 11. 

For the p-induced fission of 232 Th, the calculated ki¬ 
netic energy [closed (red) triangles] is nearly the same for 
the two studied reactions at 27 and 63 MeV. The values 
are lower than the experimental data [closed (black) tri¬ 
angles] which indicate an increase of the kinetic energy 
with the increase of the excitation energy of the fissioning 
system. For the p-induced fission of U, the calculated ki¬ 
netic energy [closed (red) circles] increases slightly when 
we go from 10 to 30 MeV. At higher energies, there is 
a small but continuous decreasing trend. The available 
experimental data |jl], [98[ for the first three proton en¬ 
ergies [closed (black) circles] are higher than our calcula- 
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tions, with an increasing trend from 30 to 60 MeV. Along 
with the proton-induced reactions studied in this work, 
we mention the recent experimental data of [113] on the 
fission kinetic energies of the neutron-induced fission of 
235 U in a broad energy range below 50 MeV. We note 
that the measured fission energies in the neutron energy 
range 30-45 MeV are approximately 162 MeV, in over¬ 
all agreement (albeit higher) with our calculation for the 
p-induced 235 U fission (Fig. 12). 

We can relate the lower kinetic energy obtained by 
CoMD, as compared to the experimental data, to the ob¬ 
servation that the CoMD code implies emission of about 
two (on average) pre-scission protons even at the lower 
energy fission reactions at 10 and 30 MeV. This is un¬ 
physical as can be concluded from our up-to-date exper¬ 
imental and theoretical understanding of the low-energy 
fission process. We discuss this feature quantitatively 
in the following. Furthermore, in all reactions studied, 
the CoMD calculations of the total energy with the soft 
symmetry potential [(blue) points connected with (blue) 
dotted lines] are lower than the corresponding ones with 
the standard symmetry potential [(red) points connected 
with (red) solid lines]. This may point to a scission con¬ 
figuration with a more elongated shape (and longer neck) 
in the soft case, as compared to a more compact shape 
(and a shorter neck) in the standard case. 


E. Pre-scission and Post-scission Particle Emission 


In Fig. 13 we show the calculated average pre-scission, 
post-scission and total proton multiplicity from the stud¬ 
ied reactions. Of course, the pre-scission proton emission 
should not be present at the lower energies, but it should 
compete with the pre-scission neutron emission at higher 
energies. A clear increasing trend is present in the calcu¬ 
lations with a substantial increrase at the highest energy 
of 660 MeV. We think that further detailed investigation 
is necessary to understand this feature of pre-scission pro¬ 
ton emission of the code at low energies (see also below 
the corresponding situation for neutrons). 

We now discuss the predictions of the CoMD code con¬ 
cerning the pre-scission, post-scission and total neutron 
multiplicity. We remind that the pre-scission neutron 
multiplicity serves as a clock of the evolution of the fis¬ 
sioning system up to the moment of scission, whereas the 
post-scission neutron multiplicity can be directly related 
to t he excitation energy of the nascent fission fragments 


114]. Both quantities can be obtained experimentally 


with proper, albeit especially difficult, measurements and 
model analysis and can offer very sensitive observables for 
any dynamical model of fission. 

From the present study with the CoMD code, we show 
in Fig. 14 the pre-scission, post-scission and total neu¬ 
tron multiplicities versus proton energy. These quantities 
show an overall increasing trend with increasing energy 
for the studied fissioning systems. 

More specifically, for the p-induced fission of 232 Th, 


the calculated pre-scission neutron multiplicity [closed 
(red) triangles] is nearly 3 at both energies 27 and 63 
MeV (showing a small increassing trend). The experi¬ 
mental value is nearly 1 at the energy of 27 MeV, and 
thus about two units lower than the calculation. This 
observation now reveals that the CoMD calculation pre¬ 
dicts a larger pre-scission neutron multiplicity (by about 
two units), as was the case for the pre-scission proton 
multiplicity that we examined before (that was assumed 
responsible for the observed lower total kinetic energy of 
the fission fragments, as compared to the experimental 
values). This conclusion, along with the corresponding 
one in regards to pre-scission proton emission, calls for 
further detailed study of the parameters of the CoMD 
code. We speculate that a careful fine-tuning of the sur¬ 
face term may be necessary to suppress the observed un¬ 
realistic feature of both pre-scission proton and neutron 
emission of the code at low energies. 

Furthermore, we note that the experimental value for 
the pre-scission neutron multiplicity at 63 MeV is in rea¬ 
sonable agreement with the CoMD calculation. In re¬ 
gards to the post-scission multiplicities, we can say that 
the calculation is in fair agreement with the available 
data, albeit larger at the higher energy. 

For the p-induced fission of 235 U, the calculated pre¬ 
scission and post-scission neutron multiplicities (Figs. 
14a and 14b, respectively) increase steadily as the en¬ 
ergy increases from 10 to 100 MeV. Agreement is seen 
with the experimental point at 63 MeV for 238 U taken 
from [2^| . For the p-induced fission of 238 U at 100 MeV, 
the pre-scission and post-scission neutron multiplicities 
are slighly higher than the corresponding values for 235 U, 
of course reflecting the larger neutron content of the for¬ 
mer nucleus. Similar observations pertain to the total 
neutron multiplicities (Fig. 14c) for the reactions stud¬ 
ied. We also note the agreement of the calculated value 
with the experimental point at 660 MeV jl05| . 


F. Fission Timescale 

We will complete the presentation of the predicted fis¬ 
sion characteristics with a discussion of the fission time 
as obtained directly by our fully dynamical CoMD cal¬ 
culations. We point out that it is a difficult task to ex¬ 
tract the fission time scale from experimental data (e.g. 
mmy Furthermore, the exctracted values are un¬ 
avoidably model and method dependent. On the other 
hand, a fully dynamical code, either a macro scop ic one, 
as the current advanced Langevin codes (e.g. [50H52j] 4 or 
a microscopic TDHF-type code (e.g. 0 ) can in princi¬ 
ple provide realistic information on the fission time scale, 
as long as, the code has been extensively benchmarked 
by comparison of its predicted fission observables with 
available experimental data on mass yield distributions, 
kinetic energies and neutron multiplicities. 

In the present study, for the first time the semiclassi- 
cal N-body CoMD code was tested with p-induced fission 






reactions and, as our presentation so far has indicated, 
it performed in an overall satisfactory manner, especially 
for intermediate and high energy fission reactions. The 
CoMD code, as a fully dynamical code, naturally de¬ 
scribes the complete dynamical path of the fission pro¬ 
cess. Therefore, with CoMD we can determine the fission 
time in a direct way. 

In Fig. 15, we show the extracted average fission time 
versus the proton energy of the reactions studied. The 
presentation of the calculations follows the same pattern 
as in Figs. 10-12. Two main groups of points are shown 
in Fig. 15. The upper group (closed points) corresponds 
to the CoMD predictions of the fission time using the full 
ensemble of the fissioning nuclei for each fission reaction. 
The lower group (open points) corresponds to the CoMD 
predictions with a selection made on the fissioning nu¬ 
cleus at the moment of fission to have exactly the initial 
Z value (thus, assuming no pre-scission proton emission). 
Within each group, the points connected with full (red) 
lines correspond to CoMD calculations with the standard 
symmetry potential, whereas the points connected with 
dotted (blue) lines are with the soft symmetry potential. 
When the selection of no pre-scission proton emission is 
made (lower group), the chosen fissioning nuclei are both 
more excited and more fissile (compared to those result¬ 
ing after the emission of pre-scission protons), and thus 
their fission time is correspondingly lower. As a gen¬ 
eral observation, within both groups, for each reaction, 
there is an overall decreasing trend with increasing pro¬ 
ton energy. Furthermore, we notice that the choice of the 
soft symmetry potential results in faster fission dynam¬ 
ics. This can be understood qualitatively by the higher 
potential energy that the soft symmetry potential implies 
for the neutron-rich low-density neck region for a highly 
deformed fissioning nucleus and, in turn, can be related 
to the corresponding lower total fission energy observed 
in Fig. 12. 

For the Th reactions, the fission time is slighly higher 
than the U reactions at nearly the same energy, reflecting 
the lower fissility of the 232 Th nucleus. We notice that 
the fission time for the p-induced fission of 235 U is on 
average longer than that of 238 U at the energy of 100 
MeV. From Fig. 13 we see that the average pre-scission 
and post-scission proton multiplicity is similar for these 
two reactions, whereas from Fig. 14 the pre-scission and 
post-scission neutron multiplicity is larger in the latter 
case. We cannot provide a simple explanation of this 
fission time difference (which is similar in both groups 
and, within each group, for each choice of the symmetry 
potential). This difference should be related to the details 
of the dynamical evolution of the two different fissioning 
systems. 

The above discussion indicates that the CoMD code 
can provide detailed information of the fission time un¬ 
der various conditions of excitation energy and fissility. 
From our study so far, if we exclude the low energy fis¬ 
sion, we expect that from 60 MeV and above, the time 
scale information may be consider realistic. It would be 


very interesting if the present predictions can be com¬ 
pared to experimental information when such informa¬ 
tion becomes available. Furthermore, it is exciting to no¬ 
tice the sensitivity of the fission time scale to the density 
dependence of the nucleon-nucleon symmetry potential, 
and thus isospin part of the nuclear equation of state, 
suggesting this observable as an additional probe of the 
latter jTlCj. 


G. CoMD Energies of Fissioning Nuclei 

After the above detailed discussion of CoMD observ¬ 
ables, we will close our presentation by examining the 
evolution of the interaction and kinetic energies of fis¬ 
sioning nuclei in their course toward scission. 

As a representative example, we show in Fig. 16 the 
results for the reaction p(30MeV)+ 235 U for the standard 
symmetry potential [(red) points connected with solid 
line] and the soft symmetry potential [(blue) points con¬ 
nected with dotted line]. In order to examine the evolu¬ 
tion of the average energies in the course to scission, and 
given the broad distribution of fission times (as seen in 
Fig. 15), in Fig. 16 the time is referenced with respect 
to the moment of scission, taken to be at t=0 fm/c. 

In Fig. 16a, the event-average nuclear interaction en¬ 
ergy is presented, taken to be the sum of the two-body, 
three-body, surface and symmetry energy terms. An 
overall increase of the interaction energy is observed as 
the fissioning system approaches the moment of scission, 
for both the standard and the soft symmetry potentials. 
Interestingly, the calculation with the soft symmetry po¬ 
tential results in higher interaction energy of the fission¬ 
ing system, as can be understood by the overall larger 
nuclear symmetry energy in the low-density neck region, 
and the overall more repulsive dymanics implied. In Fig. 
16b, the total Coulomb energy is presented, being nearly 
the same for the two choices of the symmetry potential 
(only slightly lower for the soft symmetry potential, since 
as we discussed, it results in more elongated configura¬ 
tions). A monotonic decrease of the Coulomb energy is 
observed as the fissioning system evolves toward scission, 
this decrease being essentially the main driving force of 
the nuclear system to fission. 

In Fig. 16c, the total potential energy of the fissioning 
system, namely the sum of the nuclear interaction energy 
and the Coulomb energy, is shown. In Fig. 16d, the 
kinetic energy of the fissioning system is shown, being 
lower for the soft symmetry potential, that also results 
in higher potential energy (Figs. 16a, 16c) involving more 
elongated configurations and more repulsive dynamics in 
the low-density neck region. Finally, in Fig. 16e, the total 
energy of the fissioning system is shown, being slightly 
higher for the soft symmetry potential. The decrease 
with time is due to the prescission emission of neutrons 
and protons (Figs. 13 and 14). 

In Fig. 17, we show the variation of the energies of the 
fissioning system p+ 235 U with respect to the change of 



9 


the incident proton energy from 10 MeV [(blue) points 
connected with dotted line] to 30 MeV [(red) points con¬ 
nected with solid line] to 60 MeV [(green) points con¬ 
nected with dashed line]. The calculations are with the 
standard symmetry potential. In Fig. 17a, we observe 
that with proton energies of 10 and 30 MeV, the aver¬ 
age nuclear interaction energy of the fissioning systems 
is nearly similar. The Coulomb energy (Fig. 17c) and 
the kinetic energy (Fig. 17d) increase in going from 10 
to 30 MeV proton energy. It appears that the additional 
proton energy brought in the fissioning system is stored 
as kinetic energy (Fermi motion) and Coulomb energy. 
This increase is reflected in the total energy (Fig. 17e). 

Interestingly, in going from 30 to 60 MeV proton en¬ 
ergy, we observe that the nuclear interaction energy is in¬ 
creased by nearly this amount of energy (about 30 MeV), 
whereas the Coulomb energy is nearly the same. The ki¬ 
netic energy decreases slighly, which is rather counterin¬ 
tuitive: we would expect an increase in the kinetic energy 
as we saw in going from 10 to 30 MeV proton energy. 
Thus, for the case of 60 MeV protons, the additional 
amount of energy brought in by the proton is stored 
as nuclear potential energy, indicating that essentially 
above the Fermi energy, the nuclear (mean-field) poten¬ 
tial energy is effectively momentum depend. The origin 
of this momentum dependence is in the Pauli correlations 
imposed by the CoMD procedure, i.e. the phase-space 
constraint imposed by CoMD to ensure the Fermionic 
behavior of the classically evolving system of Gaussian 
wave packets (see Section II). The total energy for the 
case of 60 MeV protons (Fig. 17e) is consistently above 
that of the previous two energies and has a diminishing 
behavior with time toward scission due to the emission 
of prescission particles, as in the other two energies. 

The above examination of the CoMD energy variations 
of the fissioning system provides a good check of the con¬ 
sistency and accuracy of the code, as applied to the de¬ 
scription of a deforming system as it evolves toward scis¬ 
sion. We wish to point out that from the present calcula¬ 
tions we cannot obtain information regarding the fission 
barriers of the involved fissioning nuclei. The reason is 
that the calculations are performed at high enough en¬ 
ergy, so that the fissioning systems are above the fission 
barrier expected to be near 6-8 MeV. In order to obtain 
average fission barriers with CoMD, a different method¬ 
ology has to be followed: the total energy of a fissioning 
nucleus has to be obtained as a function of deformation, 
placing the nucleus into a deformed harmonic potential. 
This interesting project, however, requires further com¬ 
putational effort beyond the scope of the present paper 
that we plan to undertake in the near future. 


IV. DISCUSSION AND CONCLUSIONS 

In the present work we employed the semi-classical mi¬ 
croscopic N-body code CoMD to describe proton-induced 
fission of 232 Th, 235 U and 238 U nuclei at various ener¬ 


gies. In retrospect, we chose these nuclei because of the 
availability of recent literature data and because of their 
significance in current applications of fission. We found 
that the CoMD code in its present implementation is able 
to describe fission at higher energies (i.e., above E p =60 
MeV) where shell effects are mostly washed out. We re¬ 
mind that the effective nucleon-nucleon interaction em¬ 
ployed in the code has no spin dependence, and thus the 
resulting mean field has no spin-orbit contribution. 

The total fission cross sections of the studied reactions 
were rather well reproduced. Furthermore, the ratio of 
fission cross sections over residue cross sections showed 
sensitivity to the choice of the nucleon-nucleon symme¬ 
try potential and, thus, to the density dependence of the 
nuclear symmetry energy. Consequently, this ratio can 
be used as a probe of the nuclear equation of state at low 
density and moderate excitation energy, corresponding 
to intermediate and high energy light-particle induced 
fission. Concerning total fission energies and neutron 
multiplicities, we found that they were rather adequately 
reproduced by the CoMD code (except at the lower en¬ 
ergies). Finally, information on the fission time scale can 
be obtained from the present calculations. The obtained 
fission times show a dependence on the excitation energy 
of the nucleus, as well as on the choice of the symmetry 
potential. Thus, the fission time offers one more observ¬ 
able sentitive to the isospin part of the nucleon-nucleon 
effective interaction. 

In regards to the N/Z properties of the fission frag¬ 
ments, the CoMD code appears to perform well for the 
case of p (660 MeV) + 238 U that we tested in this 
work. We note that the majority of the data of p- 
induced fission reactions have been performed in direct 
kinematics with on-line techniques and, thus, Z informa¬ 
tion of the fission fragments cannot be obtained. Apart 
from off-line gamma-ray techniques, such information has 
been obtained in stu dies in inverse kinematics mostly at 
high energies (e.g. [1071 11. In parallel to these experi¬ 

mental efforts, we mention the novel mass-spectrometric 
study of the reaction 238 U (6.5 MeV/nucleon) + 12 C in 
which proton-pickup and other channels leading t o fis - 
sion were chosen by kinematical reconstruction [ml mil. 
From an application point of view, it would be very im¬ 
portant to obtain experimental information in inverse 
kinematics at proton energies from 100 MeV to 1000 
MeV. As we mentioned earlier, we plan to perform de¬ 
tailed calculations of high-energy proton-induced fission 
in the near future. Of course the CoMD code can be 
used for neutron and other light particle incuced fis¬ 
sion reactions, heavy-ion fusion/fission reactions, as well 
as multinucleon-transfer/fission reactions and its predic¬ 
tions can be compared with existing or future experimen¬ 
tal data. We also plan to undertake calculational and 
experimental efforts in this broad direction in the future. 

To conclude, in the present study, the semiclassical 
N-body code CoMD was tested for the first time with 
p-induced fission reactions. We found that the code per¬ 
forms in an overall satisfactory manner, providing a de- 
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scription of the full dynamics of the fission process, es¬ 
pecially for intermediate and high energy fission reac¬ 
tions. We suggest that inclusion of spin dependence in 
the nucleon-nucleon effective interaction and further im¬ 
provements of the code should be implemented toward 
achieving a satisfactory description of lower energy fis¬ 
sion in which shell effects play a dominant role. We 
point out that the code parameters, as specified predom¬ 
inantly by ground-state properties of nuclei and nuclear 
matter, do not dependent on the specific reactions being 
explored and, as such, the CoMD code can offer valuable 
predictive power for the dynamics of the fission process 
in a broad range of excitation energy. Consequently, the 
CoMD code can be used for the study of fission of not 
only stable nuclei, but also of very neutron-rich (or very 
neutron-deficient) nuclei which have not been studied ex¬ 
perimentally to date and may provide guidance to up¬ 
coming RIB experiments. Moreover, this possibility can 
be further exploited in studies of fission recycling TlMltif , 
namely, the upper end of the r-process nucleosynthesis by 


the fission of the resulting very neutron-rich heavy nuclei. 
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FIG. 1: (Color online) Snapshots of the CoMD predicted 
time evolution for the reaction p (63 MeV) + 232 Th in the 
center of mass system. (Blue points: neutrons, red points: 
protons.) a) t = 0 fm/c, the proton approaches the 232 Th nu¬ 
cleus. b) t = 2000 fm/c, the fissioning nucleus is substantially 
deformed (saddle configuration), c) t = 4000 fm/c, receding 
fission fragments. For this fission event, the fission time is 
2500 fm/c. 
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p(63 MeV) + 232 Th 




FIG. 2: (Color online) CoMD calculated time evolution of the 
axial quadrupole moment Q 20 (a) and the root-mean-square 
(RMS) radius (b) of the fissisoning system from the reaction 
p (63 MeV) + 232 Th. The arrow indicates the moment of 
scission. 
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FIG. 3: (Color online) a) normalized mass distributions (cross 
sections) of fission fragments from p (27 MeV) + 232 Th. Full 
points (black): experimental data [20] . Open points: CoMD 
calculations with the standard symmetry potential, b) as 
above, but CoMD calculations with the soft symmetry po¬ 
tential. 
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FIG. 4: (Color online) a) normalized mass distributions 

(cross sections) of fission fragments from p (63 MeV) + 232 Th. 
Full points (black): experimental data [29(. Open points: 
CoMD calculations with the standard symmetry potential, 
b) as above, but CoMD calculations with the soft symmetry 
potential. 
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FIG. 5: (Color online) a) normalized mass distributions 

(cross sections) of fission fragments from p (10 MeV) + 235 U. 
Full points (black): experimental data [931 ■ Open points: 
CoMD calculations with the standard symmetry potential, 
b) as above, but CoMD calculations with the soft symmetry 
potential. 
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FIG. 6: (Color online) a) normalized mass distributions 

(cross sections) of fission fragments from p (3 0 MeV) + 235 U. 
Full points (black): experimental data [9||. Open points: 
CoMD calculations with the standard symmetry potential, 
b) as above, but CoMD calculations with the soft symmetry 
potential. 










19 



—.— * —.—.—.—i—.—.—.—i—.—.—.—i—.—.— • 

60 80 100 120 140 160 

A ff 


FIG. 7: (Color online) a) normalized mass distributions, 
cross sections of fission fragments from p (60 MeV) + U. Full 
points (black): experimental data for p (60 MeV) + 238 U [2g| . 
Open points: CoMD calculations with the standard symmetry 
potential for p (60 MeV) + 235 U. b) as above, but CoMD 
calculations with the soft symmetry potential. 
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FIG. 8: (Color online) a) normalized mass distributions 

(cross sections) of fission fragments from p (660 M eV) + 238 U. 
Full points (black): experimental data [l03f | and [104 ]. Open 
points: CoMD calculations with the standard symmetry po¬ 
tential. b) as above, but CoMD calculations with the soft 
symmetry potential. 
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FIG. 9: (Color online) a) Fission fragment average Z vs A 
from p (660 MeV) + 238 U. Full points (black): experimental 
data Il03l . 1 1041 . Solid (red) line: CoMD calculations with the 
standard symmetry potential. Dotted (blue) line : CoMD cal¬ 
culations with the additional selection of the fissioning system 
not to emit any pre-scission protons (Z=93, see text), b) stan¬ 
dard deviation of the isobaric Z distribution vs A. Symbols 
and lines as above. 
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FIG. 10: (Color online) Calculated total fission cross section 
with respect to proton energy E p . The CoMD calculations 
are carried out with the standard symmetry potential and 
are shown with full (red) symbols connected with full (red) 
lines. The reactions are indicated as follows: triangles: p (27, 
63 MeV) + 232 Th, circles: p (10, 30, 60, 100 MeV) + 235 U, 
squares: p (100, 660 MeV) + 238 U. Some experimental data 
are shown with closed (black) symbols as follows: triangles: 
p(27, 63 MeV) + 232 Th gjj, circles: p (63 MeV) + 238 U 
[91J |. square: p (660 MeV) + 238 U [103 . 1041 ]. The point at 
E p =660 MeV is displayed at E p =160 MeV. The errorbars are 
statistical. 
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FIG. 11: (Color online) Calculated ratio of the fission cross 
section over residue cross section with respect to proton en¬ 
ergy. CoMD calculations with the standard symmetry po¬ 
tential are with full (red) symbols connected with full (red) 
lines. Calculations with the soft symmetry potential are with 
full (blue) symbols connected with dotted (blue) lines. The 
reactions are indicated as follows: triangles: p (27, 63 MeV) 
+ 232 Th, circles: p (10, 30, 60, 100 MeV) + 235 U, squares: 
p (100, 660 MeV) + 238 U. The points at E p =660 MeV are 
displayed at E p = 160 MeV. The errorbars are statistical. 
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FIG. 12: (Color online) Calculated average total energy of 
fission fragments with respect to proton energy. CoMD cal¬ 
culations with the standard symmetry potential are with full 
(red) symbols connected with full (red) lines. Calculations 
with the soft symmetry potential are with full (blue) symbols 
connected with dotted (blue) lines. The reactions are indi¬ 
cated as follows: triangles: p (27, 63 MeV) + 232 Th, circles: 
p (10, 30, 60, 100 MeV) + 235 U, squares: p (100, 660 MeV) + 
238 U. Some experimental data are shown with closed (black) 
symbols as follows: triangles: p (27, 63 MeV) + 232 Th [29ll . 
circles: p (10, 30 MeV) + 235 U [jj| and p (63 MeV) + 238 U 
[9l| . The points at E p =660 MeV are displayed at E p =160 
MeV. (Statistical errorbars on the theoretical points are about 
1.0-1.5 MeV and are not shown for clarity.) 
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FIG. 13: (Color online) Calculated average proton multiplic¬ 
ities with respect to incident proton energy: a) pre-scission b) 
post-scission and c) total multiplicities. CoMD calculations 
are carried with the standard symmetry potential. The reac¬ 
tions are indicated as follows: triangles: p (27, 63 MeV) + 
232 Th, circles: p (10, 30, 60, 100 MeV) + 235 U, squares: p 
(100, 660 MeV) + 238 U. The point at E p =660 MeV is dis¬ 
played at E p =160 MeV. 
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FIG. 14: (Color online) Calculated average neutron multi¬ 
plicities with respect to incident proton energy: a) pre-scission 
b) post-scission and c) total multiplicities. CoMD calculations 
are carried with the standard symmetry potential. The reac¬ 
tions are indicated as follows: triangles: p (27, 63 MeV) + 
232 Th, circles: p (10, 30, 60, 100 MeV) + 235 U, squares: p 
(100, 660 MeV) + 238 U. Some experimental data are shown 
with closed (black) symbols as follows: triangles: p (27, 63 
MeV) + 232 Th [2S|, circle: p (63 MeV) + 238 U [2<|. square: 
p (660 MeV) + 238 U [HJ. The point at E p =660 MeV is 
displayed at E p =160 MeV. 
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FIG. 15: (Color online) Calculated fission time with respect 
to incident proton energy. CoMD calculations with the stan¬ 
dard symmetry potential are with (red) symbols connected 
with full (red) lilnes. Calculations with the soft symmetry po¬ 
tential are with (blue) symbols connected with dotted (blue) 
lines. The full symbols (upper half of the figure) are with 
the full ensemble of the fissioning nuclei, whereas the open 
symbols (lower half) are with the selection of the fissioning 
system not to emit any pre-scission protons (see text). The 
reactions are indicated as follows: triangles: p (27, 63 MeV) 
+ 232 Th, circles: p (10, 30, 60, 100 MeV) + 235 U, squares: 
p (100, 660 MeV) + 238 U. The points at E p =660 MeV are 
displayed at E p = 160 MeV. 
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FIG. 16: (Color online) Evolution of ensemble average CoMD 
energies of fissioning nuclei in the interaction of p(30 MeV) 
with 235 U. The moment of scission is taken as t=0 fm/c. 
(Red) points connected with solid lines are with the stan¬ 
dard symmetry potential. (Blue) points connected with dot¬ 
ted lines are with the soft symmetry potential. The energies 
are: a) Nuclear potential energy (the sum of the two-body 
interaction energy, the three-body interaction energy, the sur¬ 
face energy and the symmetry energy), b) Coulomb potential 
energy, c) Total potential energy [sum of energies of (a) and 
(b)[. d) Kinetic energy, e) Total energy [sum of (c) and (d)[. 
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FIG. 17: (Color online) Evolution of ensemble average CoMD 
energies of fissioning nuclei in the interaction of p with 235 U. 
The reactions are with: 10 MeV protons: (blue) points con¬ 
nected with dotted line, 30 MeV protons: (red) points con¬ 
nected with solid line, 60 MeV protons: (green) points con¬ 
nected with dashed line. The moment of scission is taken as 
t=0 fm/c. The calculations are with the standard symmetry 
potential. The energies are: a) Nuclear potential energy (the 
sum of the two-body interaction energy, the three-body inter¬ 
action energy, the surface energy and the symmetry energy), 
b) Coulomb potential energy, c) Total potential energy [sum 
of energies of (a) and (b)]. d) Kinetic energy, e) Total energy 
[sum of (c) and (d)]. 







